Suppression of exciton dephasing in quantum dots through ultrafast multipulse control 
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We investigate the usefulness and viability of the scheme developed by Viola and Lloyd [Phys. Rev. A 58, 
2733 (1998)] to control dephasing in the context of exciton-based quantum computation with self-assembled 
quantum dots. We demonstrate that optical coherence of a confined exciton qubit exposed to phonon-induced 
dephasing can be substantially enhanced through the application of a simple periodic sequence of control pulses. 
The shape of the quantum dot has a significant effect on the dephasing properties. Remarkably, we find that 
quantum dots with parameters optimized for implementing quantum computation are among the most suscep- 
tible to dephasing, yet periodic decoupling is most efficient for exactly that type of dot. We also show that 
the presence of an electric field, which is a necessary ingredient for many exciton-based quantum computing 
schemes, may further increase the control efficiency. Our results suggest that dynamical decoupling may be a 
method of choice for robust storage of exciton qubits during idle stages of quantum algorithms. 

PACS numbers: 03.67.Lx,73.21.La,81.07.Ta 



c 

o 
o 



> 

cn 
o 

oo 
O 



X 



I. INTRODUCTION 

The quest for practical implementations of quantum infor- 
mation processing (QIP) 1 has resulted in a wide variety of pro- 
posed quantum computing architectures, based on virtually 
all kinds of different physical systems 2 -. Thanks to continu- 
ous advances in nanoscale design and manufacturing of low- 
dimensional structures, solid-state platforms promise, in par- 
ticular, a high level of scalability and large-scale integration. 
Among the various proposals which exist to date, schemes 
which rely on charge degrees of freedom of carriers (exci- 
tons) confined in a self-assembled semiconductor quantum 
dot (QD) 3 are especially appealing in view of current ultrafast 
spectroscopic capabilities, which make it possible to access 
sub-picosecond gating times by means of suitable all-optical 
control techniques 4 ^. 

In practice, such advantages are hindered by the rate at 
which quantum information stored in exciton qubits and/or 
exciton ancillary states is irreversibly lost due to dephasing 
- as characterized by a typical time scale T^. For QDs to 
be employed in quantum information devices, decoherence 
processes, which involve energy exchange with the host crys- 
tal are highly suppressed due to the large splitting between 
quantized energy levels. Thus, the dominant contribution to 
Ti arises from energy-conserving transitions which do not 
change the occupations of the logical states, but result in 
pure dephasing and unrecoverable loss of phase information. 
While typical T% times i.e. finite qubit lifetime via exciton 
recombination, for single QDs are in the nanosecond range, 
it has been shown theoretically 7 ' 8,9 and experimentallyiS that 
substantial dephasing can occur in much shorter time scales 
(one or two picoseconds at temperatures of a few K), mainly 
due to coupling with acoustic phonons. 

Borrowing inspiration from coherent averaging tech- 
niques in high-resolution nuclear magnetic resonance 
spectroscopy 11 , dynamical decoupling (DD) methods have re- 
cently emerged in QIP as a versatile strategy for counter- 
acting decoherence from non-Markovian environments, in a 
way which is substantially less resource-demanding than tra- 



ditional quantum error-correcting codes. In its simplest form, 
DD is implemented by subjecting the system to periodic se- 
quences of fast and strong (so-called bang-bang) control op- 
erations drawn from a basic repertoire, so that interactions 
whose correlation times are long compared to the control time 
scale are effectively averaged out to finite accurac y 12 ' 13 . In 
particular, a number of analytical and numerical studies have 
established the potential for DD schemes based on high-level 
concatenated 14 or randomized 1516 design to extend the coher- 
ence time of spin degrees of freedom in gate-defined semi- 
conductor QDs by (at least) two orders of magnitud e) 17 ' 1 8 i 1 9 i 20 . 
Remarkably, simple DD protocols consisting of a single or 
multiple spin-echoes have been experimentally demonstrated 
in (gate-defined) double QD devices 21 and rare-earth solid- 
state centers 22 , respectively. In this paper we will consider 
self-assembled semiconductor QDs. 

The usefulness of a suitable burst of ultrashort pulses as a 
tool to reduce the initial optical polarization dephasing in this 
type of QD has been demonstrated by Axt et a££. In this work, 
we investigate the effectiveness of multipulse DD toward both 
mitigating the initial coherence decay and ensuring enhanced 
coherence preservation over a desired storage time scale, with 
emphasis on QIP applications. Our analysis has two main im- 
plications: first, we show that in spite of being most prone 
to dephasing, QDs which are optimized to support quantum 
computation are also most efficiently stabilized by DD. Fur- 
thermore, thanks to the fact that the exciton spectral density is 
down-shifted in frequency by the application of a static elec- 
tric field, DD performance are further improved for a biased 
QD as required for most computational applications 4 ^. 

The content of the paper is organized as follows. The next 
section provides an in-depth analysis of the free dephasing dy- 
namics of an exciton qubit in the absence of applied control. 
After laying out the relevant QD model in Sec. IIA, special 
emphasis is devoted to isolate and analyze the contributions 
of the piezoelectric and deformation couplings to the exciton 
spectral density (Sec. IIB) for various dot shapes, as well as 
to quantify the influence of the spectral density details on the 
dephasing behavior. In Sec. Ill, we summarize our results on 
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the controlled dephasing dynamics, and characterize the de- 
pendence of DD performance upon different system and con- 
trol parameters - including dot shape, temperature, separation 
between pulses, and application of an external bias field. In 
particular, DD performance for robust exciton storage in QIP 
applications are highlighted in Sec. IIIC, HID, respectively. 
Concluding remarks follow in Sec. IV. 



H. DEPHASING DYNAMICS OF CONFINED EXCITONS 
A. QD model Hamiltonian 

We consider a GaAs/AlAs QD with either or 1 ground 
state exciton 9 , which corresponds to the qubit logical states 
|0) or |1) respectively. The Hamiltonian of this two-level ex- 
citonic system interacting with the phonon modes of the lattice 
is given by 

H = E exc c^c+hY^^(^)b} M b jM +h^cJ2(9ikbl k +H.c.), 

(1) 

where E exc is the energy of the ground state exciton relative 
to the crystal ground state, c^(c) are creation (annihilation) 
operators for an exciton, frt k (6j k) are bosonic creation (anni- 
hilation) operators for a phonon of mode j, wave vector k, and 
angular frequency ujj (k), and gj^ is the coupling between the 
exciton and a phonon of mode j, k, respectively. 

We assume that the exciton and the phonon bath are initially 
uncorrelated, with the phonon bosonic reservoir being in ther- 
mal equilibrium at temperature T. As time evolves, the qubit 
becomes entangled with the environment and the off-diagonal 
elements of the exciton density matrix evaluated at time t in 
the interaction picture with respect to the free system and bath 
Hamiltonians are given by^i 2 - 



be safely neglected. For acoustic phonons, we assume the 
dispersion relation Wj(k) = Wj|k|. In addition, we assume 
a small lattice mismatch at the boundary between the QD 
and the barrier, and therefore approximate the coupling to the 
phonon bath as the bulk coupling modulated by the appropri- 
ate form factor—: 

<&, k = j d 3 r e d 3 rfc|*(r e ,r h )| 2 (G? k e ik-re - Gj k e ik - rh ). 

(5) 

Here, ^(r e , r/j) is the exciton wave function and Gj£ is the 
bulk coupling of the single particle to the phonon bath, 
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v / 2ghu J (k)V 



[kD e / h + iM 3 (k)}. (6) 



The index j runs over the two transverse and the longitudinal 
modes, g is the density of the QD, and V is a normalization 
volume. Excitons in bulk semiconductor couple to longitu- 
dinal acoustic phonons through a deformation potential cou- 
pling D e J h , and to all the modes considered through a piezo- 
electric potential coupling Mj . For zincblende crystals such 
as GaAs/AlAs, 



M,(k) 



2ei4e 



(kkytP + + kki^)- (7) 



In what follows, we shall approximate \Mj\ 2 by its angular 



average 
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d0sin(0)M 2 (k) = Aj 
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(8) 



Poi(t) = P* 10 (t) = poi(0)e 



-r(t) 



(2) 



where e is the electron charge, is the piezoelectric coeffi- 
cient, and e s , eq denote the relative permittivity of the material 
and of free space, respectively. £ x / y / z are unit vectors describ- 
ing the polarization of mode k, and Aj are mode-dependent 
geometrical factors. 



with 



r(t) 



, i(uj) , / hut 

dui — 7t- cotn 



(^)(l-co sM) ). (3) 



Here, ks is the Boltzmann constant, and the spectral function 



8(w -^(k))!^,!, 



(4) 



describes the interaction of the qubit with a phonon of fre- 
quency ojjQt), For an exciton qubit, T(t) is formally a factor 
of 4 smaller than for a spin coupled to a bath of harmonic os- 
cillators as described in Ref. [12J. This is due to the different 
Hamiltonians. For the exciton qubit the coupling to the envi- 
ronment only occurs when the qubit is in the logical state |1) 
(exciton present), while a spin qubit couples to the environ- 
ment in both logical states. 

In order to maximize coherence, we consider sufficiently 
low temperatures, so that the coupling to optical phonons may 



We will consider nanostructures in the strong confinement 
regime, so that Coulomb coupling may be neglected and ex- 
citons may be accurately modeled as products of Gaussian 
wave-packets within the single particle approximation. We 
assume that the QD is isotropic in the plane perpendicular to 
the growth direction. The excitonic wave-function is then 



(9) 



where 



/^ 2 ) : 



■ exp 



e/h 



-e/h 



2\l 



2A? 



(10) 

and2A z .. and2A r ,, are the characteristic widths of the elec- 
tron/hole wavefunction in the z (growth) direction, and in the 
plane r = ^/ x 2 + y 2 respectively, and are related in the usual 
way to the characteristic frequencies of the parabolic confin- 
ing potential in the appropriate direction. For later reference, 
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we define a characteristic volume of the QD, 



V, h 



8A A 2e 



(11) 



For the exciton-phonon system, Eq. can be written as 

I exc {uj) = I e {uj) + I h {uj) - 2I eh (u>), (12) 

where I e /h (ui) is the spectral density of the electron/hole, and 
I e h (w) modifies the spectral density due to the electron hole 
interference. Each term has the following form: 



/ S H~5>|G;->)| 2 / S (^)exp( 

j 3 

where the index s runs over e, h, and eh, 



x 2 y 

2v 2 



(13) 



and 



\2 

\r/z) Bh 



\ 2 -J- \ 2 

A (r/z)e + \r/ Z ) h 



(14) 



(15) 



The form of f s (x) depends on the dot shape. For X rs > X z 
(oblate), 



/f>) 



i eif(ix) 



(16) 



where erf (x) is the error function; for A rs < X Zs (prolate), 

erf (x) 



zr°(x) = 

and for A rs = X Zs (spherical) 

ft ph {x) 



(17) 



(18) 



In Eq. ( Tot . \G'j s {oj)\ 2 includes terms due to both the defor- 
mation potential and the piezoelectric coupling, that is, 



l^»| 2 = 72 ^ 



\M. 



(19) 



where (D eh ) 2 = D e D h . This allows the spectral density to 
be separated into the sum of the spectral densities related to 
each of the two coupling mechanisms, 



7(w) = I def {uj) + I piez (uj), 



(20) 



and the total evolution in Eq. (O to be expressed as a product 
of contributions from the deformation potential and piezoelec- 
tric interactions, 



e -r(t) = e -r<"f(t) e -r**"(t)_ 



(21) 



3 



3 




FIG. 1: Piezoelectric Ig iez (tj) (top) and deformation (u) (bot- 
tom) contributions to the electron spectral density for the cases of: a) 
spherical QD (A rc = 3.84 nm), b) oblate QD (A rc = 6 nm), and c) 
prolate QD (X r<s — 2 nm). Inset: In I (lu) for a single electron with 
Xr e =3.84 nm (solid line). The dashed line is the parabolic fit with 
the curve (29.7 - 2.79x 10" 25 w 2 ), see Eq. g3). 



B. Single particle and excitonic spectral density 

The dephasing dynamics of the exciton qubit is determined 
by the exciton spectral density. In general, I(u>) is positive, 
goes to zero as ui — > 0, and becomes negligible for large ui. 
Aside the material parameters, the main factor influencing the 
spectral density is the shape of the QD 7 , which enters our cal- 
culations through the values of Xr r i z \ e and A( r / Z ) h . 

To understand the system at hand, we will first consider 
the behavior of the spectral density for a single trapped par- 
ticle. Throughout this paper we will consider QDs with 
V^ e = 455.28 nm 3 . For a single particle, for instance an 
electron, the maximum of the spectral density is greatest for a 
spherical QD. For oblate or prolate QDs with the same char- 
acteristic volume V[f,, the spectral density has a lower max- 
imum, but extends to higher frequencies. The top panel of 
Fig. [1] shows the piezoelectric spectral densities for a single 
electron, If lez , for different dot shapes. As the dot becomes 
more asymmetric, the interplay between f e (oj^ e /vj) and the 
Gaussian term in Eq. ( fl~3"l ) reduces the maximum of /P ie2 (a;), 
whilst increasing I p%ez (lS) at high frequencies. 

For the deformation potential coupling, the single particle 
spectral densities follow a similar trend, however the extra lo 2 
term in the expression for lf e * (u>) (see Eq. ( fT9l )) results in 
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larger lf e ' (oj) at high frequencies. The bottom panel of Fig. Q] 
shows the deformation potential spectral densities for a single 
electron in QDs with different shapes. 

We can define an upper cut-off frequency ui e , such that I (us) 
is negligible for u> 3> oj c . The QD spectral density is com- 
posed of a sum of terms, each containing a Gaussian function. 
The high frequency behaviour of the overall function is there- 
fore very similar to a Gaussian, thus we may write 



I(u) « F{uj)e 



(22) 



and treat lo c as the cut-off frequency. In order to obtain a value 
for u> r , we consider 



ln/(w)=lnFH-^, 



(23) 



and fit this form to our numerical results. The inset of Fig. 03 
shows In I(u>) for a single electron in a spherical QD (A re = 
3.84 nm), along with fitted curve 29.7 - 2.79x10~ 25 cj 2 which 
leads to ui c = 1.89xl0 12 rads" 1 (Tiuj c = 1.25 meV). 

For an exciton in a self-assembled QD, the electron and 
hole confinement in the z (growth) direction is due to the 
band-offset between the barrier and dot layers. Because it 
is then mainly determined by the thickness of these layers, 
A 2e ph X Zh ■ The radial in-plane confinement, however, may 
be different for electron and hole, and is accurately modeled 
by harmonic potentials. This leads to 



^r h — 



i m*w e 
m* h uj h 



(24) 



where ui e / h is the characteristic frequency of the confining 
potential for the electron/hole, and m* , h is the electron/hole 
effective mass. Throughout this paper, we shall discuss GaAs 
QDs 6 , the material parameters of which are listed in Table Q] 
Because of Eq. d24l >. the electron and hole wave-functions are 
not simultaneously spherical for a given QD. There are three 
possible regimes for the exciton wave-function. For a heavily 
prolate QD, both the electron and hole wave-functions are pro- 
late (prolate-prolate regime), similarly for heavily oblate QDs 
the wave-functions are both oblate (oblate-oblate regime). In 
between these two situations, the electron and hole will stretch 
in different directions (oblate-prolate regime). The most sym- 
metric configuration in the oblate-prolate regime occurs when 
£e = £h> corresponding to the following radial confinement: 



4 A /2(1 



(25) 



In the results plotted throughout the paper, we shall assume 
TiL0 e /hu!h = 1.25, which corresponds to the parameters for 
the optimal dot for quantum computation discussed in Ref. 

For the excitonic system, there is an additional difference 
between the spectral density components due to piezoelectric 
and deformation potential couplings: for piezoelectric cou- 
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0.34 


Longitudinal geometrical factor 


A L 


3/35 


Transverse geometrical factor 


A T 


1/21 



TABLE I: Material parameters for GaAs QDs, from Ref. mo 
denotes the free electron mass. 



pling, the interference term is positive, whereas I fj is 
negative. For deformation potential coupling, the interference 
term is of the same order of l£ e f and if^^ . As a consequence, 
^exc behaves in a similar way to the single particle spectral 
density in that it has a lower maximum in the prolate-prolate 
and oblate-oblate regimes compared to QDs with A r<! » A,, e , 
however 7™£ will be greater for large u>. 

Fig. |2] shows the deformation potential spectral density 
lexl f° r different dot shapes. Notice that at A rc = A rc , the 
electron wave-function is as oblate as the hole wave-function 
is prolate; for deformation potential coupling, I]~f behaves 
similarly for A, e > A re as for A re < A rc . 

For piezoelectric coupling, the different sign of the interfer- 
ence term 1 V V^ Z leads to a very different behavior of the spec- 



tral density. First, if the electron and hole wave-functions were 
exactly the same, there would be no piezoelectric coupling, 
that is, I&f = 0. At low frequencies, Ii iez (uj) w If t iez (uj), 
and their sum cancels with the contribution from the interfer- 
ence term, leading to If^ ~ 0. The width of this deple- 
tion of the spectral density is greater for smaller A re . For 
dots in the oblate-oblate regime (A rc > 4.85 nm), as A r - e 
is decreased the maximum of the piezoelectric exciton spec- 
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FIG. 2: Excitonic deformation spectral density 1^1 (w) for the cases 
of: a) X r> . = A rc = 4.16 nm, b) oblate-oblate regime (A rc = 6 nm), 
and c) prolate-prolate regime (A rc = 3 nm). 
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Oblate -prolate and 
prolate -prolate regimes 




2 3 4 

Bco (meV) 

FIG. 3: Top panel: Excitonic piezoelectric spectral density I^ e c z (u>) 
in the oblate-oblate regime, for the cases of increasing A re . a) A,- s = 
6 nm, b) A re = 9 nm, c) A rc = 13 nm, and d) A r< . = 18 nm. Bottom 
panel: I^.% z (u)) in the oblate-prolate and prolate-prolate regimes, 
respectively, for decreasing A r<! . a) A r< , = 4 nm, b) X Tc = 3 nm, c) 
A re = 2 nm, and d) A re = 1 nm. 



tral density increases, however it also shifts towards higher 
frequency modes as shown in the top panel of Fig. [3] As 
A re decreases further, we first enter the oblate-prolate regime 
(3.84 nm< A rc < 4.85 nm) and then the prolate-prolate 
regime (A re < 3.84 nm). As A rc reduces to even smaller 
values, the piezoelectric spectral density becomes small for 
all u>. This is due to the type of confinement in the growth 
direction, which results in A 2e = X Zh . In this limit, £ e w £/j 
and If lez + I^ lez + I^ ez w for all 10. This behavior may 
be seen in the bottom panel of Fig. [3] in which If!Jf z is plot- 
ted for values of decreasing A rc . We shall denote the value of 
A re which yields the largest maximum of Tg£? as \ p r lez , with 
~KT Z Ar e . For V^ e = 455.28 nm 3 , \f e ez = 5.58 nm. The 
discrepancy between A re and \ P l lez is due to the constraint 
Az. = A z , . 



Again, we may follow the same method as in the single- 
electron case to define and estimate a spectral cut-off uj c . 
For an exciton with = A re = 4.16 nm, this yields 
u> c = 3.52xl0 12 rads -1 (hu> c — 2.32 meV), i.e. lo c is of 
the same order as the single particle case for a spherical QD. 
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FIG. 4: Decoherence dynamics for A re = 6.16 nm at 4K (solid line), 
and asymptotic decoherence value from Eq. J26b (dashed line). Inset: 
Decoherence dynamics for A rc = 6.16 nm at 77K, using: a) exact 
results, and b) the short-time approximation from Eq. d27b . 



C. Effect of spectral density on exciton dephasing 

Insight into the dephasing properties may be obtained by 
examining the evolution of the off-diagonal density matrix 
elements with time, Eq. (0. Physically, the coherence ele- 
ment poi(i) is proportional to the system optical polarization 
at time t, P(i) (see e.g. Eq. (9) in Ref. 7). In order to em- 
phasize this relation, we have chosen to monitor the quan- 
tity | e - r( *)| 2 = \P{t)\ 2 /\P(t = 0)| 2 . Fig. shows the de- 
coherence factor | exp(— r(t))| 2 for A,. e = 6.16 nm at 4K. 
| exp(— r(£))| 2 falls quickly from unity before saturating at a 
temperature- and shape- dependent value. For the examples 
in Fig. ID | exp(— r(oo))| 2 is non zero, but at higher tempera- 
tures and for differently shaped dots poi ma Y fall t° zeroi 

More specifically, in the limit of t — > oo, Eq. (f3]) becomes 



r(oo) 



. J(w) . 
aid — — cotn 



(i 



(26) 



therefore the non-diagonal density matrix may saturate to a 

finite value if I(u>) lu 3 . This was discussed for spe- 

cific spectral density functions in Refs. M and 1241 Saturation 
at a non-zero value does not occur for the spin-boson sys- 
tem with Ohmic spectral density analyzed in Ref. [13. Fig. 
[5] shows the asymptotic coherence value | exp(— r(oo))| 2 as 
a function of A re for an exciton at 77 K in the presence of 
a) deformation potential coupling alone, b) piezoelectric cou- 
pling alone, and c) both effects included. For the piezoelec- 
tric interaction, r piez (oo) — > as A, e — > due to the fact 
that If™ -> (Fig. |3j. As A re increases, TP iez (oo) be- 
comes large, reflecting the fact that the asymmetry between 
electron and hole wave-functions becomes relevant, therefore 
increasing the piezoelectric coupling. For the deformation po- 
tential coupling, r de ^(oo) is maximum very close to A re . It 
can be seen in Fig. [5] that in order to minimize the long- 
term decoherence, we should either consider very prolate QDs 
(A re <§C A 2c ) - for which piezoelectric coupling can be ne- 
glected - or consider the local maximum 

of | e -r(oo)|2_ which 
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X r (nm) 



FIG. 5: Asymptotic decoherence factor | exp(— r(oo))| 2 versus \ r< , 
for: a) deformation potential coupling, b) piezoelectric coupling, and 
c) both effects included, for an exciton at T = 77K. 
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FIG. 6: Dependence of fixed-time dephasing due to piezoelectric 
coupling, | exp(— T plez (t = t'))\ 2 , with respect to A re for an ex- 
citon at 77K. a) t' = 0.3 ps, b) t' = 3 ps, c) t' = 6 ps. In d), the 
long-time decay | exp(— r plez (oo))j 2 is also depicted. 



for low temperatures the fastest dephasing (that is, the min- 
imum width of the Gaussian in Eq. d27| >) occurs for A rc w 
A^* ez : at 4K, for instance, the fastest dephasing happens at 
X ra « 6.18 nm. For more oblate/prolate QDs, the initial de- 
phasing occurs slower. Fig. |6]shows | exp(— T piez (t = t'))\ 2 
evaluated at different evolution times as a function of A re 
and at T=77K. The short-time dephasing (t' < 0.3 ps) be- 
comes negligible for small A rc as I plez — > 0, and also be- 
comes negligible as A rc — > oo. Depending on A re , de- 
phasing may be not mono tonic with time. For example, in 
Fig. |6]exp(-rP iez (6 ps)) < exp(-r^ ez (oo)) for 10 nm 
^ A re < 33 nm. This is due to the existence of a local 
minimum in the evolution of exp(— T piez (t)) before satura- 
tion occurs 7 (see Fig. |4|. 

For deformation potential coupling, the width of the Gaus- 
sian in Eq. d27b is locally maximum for QDs with A re close 
to A rc . At 4K, this occurs at A re « 4.3 nm. The short-time 
dephasing becomes negligible for very large and very small 
X rc (see Fig. UJ where t = 0.3 ps). For low temperatures, a 
local minimum of the dephasing occurs near to A ro . As the 
temperature is increased, the local minimum is flattened and 
becomes a global maximum (see Fig. [7]). 
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occurs at r piez s» T de f , at A,. e ps 25 nm in the case of Fig. [5] 
The results in Fig. [5] are calculated at a temperature of 77K. 
The magnitude of the dephasing can be reduced by consider- 
ing lower temperatures. 

Let us now consider the opposite time limit. The short- 
time decoherence behavior offers insight into how the spec- 
tral density affects the initial dephasing rate, that is, how fast 
exp(— r(i)) falls from unity. For small t, Eq. (0 becomes 



FIG. 7: Dependence of short-time dephasing due to deformation cou- 
pling, | exp(— F de f (t — 0.3ps))| 2 , upon A re for an exciton at differ- 
ent temperatures: a) 4K, b) 15K, c) 30K, d) 56K, and e) 77K. 



(27) 



Eq. d27l > shows that at short times exp(— T(t)) has a Gaussian 
shape with a width that depends on the structure and mate- 
rial parameters, as well as temperature. The inset of Fig. [4] 
compares the approximation of Eq. (l27l i with the the exact 
dephasing curve, Eq. (01. 

For the piezoelectric interaction, the rate at which 
exp(— T piez (t)) falls from unity is temperature-dependent but 



III. PERIODICALLY CONTROLLED EXCITON 
DEPHASING DYNAMICS 



Having clarified the essential features of the phonon- 
induced dephasing dynamics in the absence of control, we 
now show how it can be substantially reduced through DD. 
Our main goal here is to assess the benefits resulting from the 
application of DD in its simplest periodic form, involving re- 
peated bit-flips as in Ref . 1 1 2l 
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A. Control setting 

Let DD be implemented by subjecting the exciton to a 
train of uniformly spaced ideal 7r-pulses, applied at instants 
ti = I At, £ = 1,2,..., with At > being the separation 
between consecutive pulses. T c = 2At defines a complete 
DD cycle, which brings the control propagator back to unity. 
The controlled dephasing dynamics may then be described by 
a modified decoherence function^, 

r(i) = ^^) coa (^)„ ()) ^(^), 

(28) 

where time is understood to be stroboscopically sampled at in- 
teger multiples of the control cycle time, t = NT C = 2N At. 

In practice, the required control rotations may be effected 
through suitable laser pulses. On resonance, the exciton qubit 
will Rabi-oscillate between the logical states, and a pulse of 
appropriate length will perform a bit-flip. Eq. d28l l assumes 
that such bit-flips are instantaneous, which is an adequate ap- 
proximation provided the time necessary to perform each bit- 
flip operation is much shorter than any other time scale rele- 
vant to the problem. In modeling the evolution of the excitonic 
qubit, however, a constraint must be placed on the timing of 
the control pulses. The central energy of a control pulse must 
be resonant with the energy of the ground state exciton, and its 
frequency spread must be small enough such that the probabil- 
ity of exciting higher energy levels is negligible. This places 
a lower limit on the pulse length, as too short a pulse would 
excite higher energy levels of the QD, and in extreme circum- 
stances ionize the electron. To account for this in a model 
which assumes instantaneous bit-flips, the delay between bit- 
flips must be significantly longer than the minimum length of 
each control pulse. 

The constraint on the pulse length is also a requirement in 
QD-based quantum computation schemes which are based on 
the so-called bi-excitonic shift 6 to perform two-qubit gates. 
The bi-excitonic shift can be of the order of a few meV 6 - , thus 
the constraint translates into a pulse length of the order of a 
fraction of a picosecond. The control pulses we are consid- 
ering in this paper, however, can be shorter. The relevant en- 
ergy scale in this case is the difference between the ground 
state and higher energy levels, and in the strong confinement 
regime of interest here this can be of the order of tens of meV. 
This translates to a minimum pulse length of the order of a few 
tens of femtoseconds, and a constraint on the interval between 
control pulses of the order of a few tenths of a picosecond. 
Finally, we assume that the lattice relaxation time is much 
greater than the interval between pulses, and that no system- 
atic and/or random control errors affect the DD operations. 



B. Dephasing suppression through multipulse control 

Figure [8] shows the decoherence factor, | exp(— T(t))\ 2 , for 
A re = 6.16 nm at 77K in the presence of a DD sequence 
with different pulse delays, whereas the inset compares the 
free evolution with the controlled evolution on a smaller time 



l r a i i i i I 
i 

0.8 - 

b 

\ ' - , t 

M _ 0.6 -\ 

£" \ M _ 1 \. J... I ...'...! . j — . 

- 4 \ C 0.6 \ 

C \ [o 0.4 \ 

o.2- x °j \- \ ■ ■ ■ i 

^, u 2 4 6 8 10 
t (ps) 

o 1 1 1 — f ---- 1 

50 100 150 200 250 

t (ps) 

FIG. 8: Decoherence dynamics, | exp(— r(t))| 2 , for A re = 6.16 nm 
at 77K under periodic DD with: a) At = 0.2 ps, b) At = 0.25 ps, 
and c) At = 0.3 ps. Inset: free (solid line) and controlled (dotted 
line) evolution for the case of At = 0.25 ps. 



scale for the case of At = 0.25 ps. Decreasing the pulse sepa- 
ration At decreases the dephasing for two reasons: firstly, the 
qubit dephases slower in the presence of more closely sepa- 
rated control pulses, and secondly the control sequence begins 
at an earlier instant, thus less coherence is lost before the first 
control pulse occurs. 

As expected, the dephasing follows the free evolution un- 
til the first bit-flip occurs, after which a transient regime sets 
in during the first few control cycles. Eventually, the con- 
trolled coherence decays away to zero, but over a substan- 
tially longer time scale than the one associated with the free 
evolution. Our numerical results suggest that no saturation at 
a non-zero value takes place for the controlled evolution in the 
relevant parameter regime, however the coherence decay time 
increases with decreasing At. This is apparent when com- 
paring the controlled evolution over different time scales, e.g. 
compare the dotted lines in the main panel and the inset of Fig. 
[S] respectively, which show the evolution for At = 0.25 ps. 
As long as DD is turned on early enough, there exists a time 
interval in which the controlled evolution dephases less than 
the free evolution. This time interval can be orders of magni- 
tude greater than the free evolution dephasing time. 

Beside depending on the control time scale At, the DD ef- 
ficiency is influenced by the shape of the dot, which enters 
Eq. (l28T > through the spectral density I(u>) analyzed in Sec. 
IIC. Alternatively, one may regard the controlled decoherence 
behavior as being determined by a DD-renormalized spectral 
density I(ui) tan 2 {u> At/2). The last term presents a reso- 
nance at Lo res = TV / At, which is the characteristic frequency 
introduced by the periodic pulsing. The smaller At, the higher 
this characteristic frequency is. Although the time-dependent 
contribution prevents the integrand from diverging to infin- 
ity, the location of cu res relative to the bare spectral density 
plays a crucial role in determining the dephasing properties: 
at sufficiently long time t, the largest contribution to T(t) for 
our system originates from frequencies close to ui res . There- 
fore, DD is expected to be efficient provided that uo res > u> c , 
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FIG. 9: Deformation decoherence factor | exp(— Y de ^ (t'))\ 2 for 
an exciton at 4K (top panel) and piezoelectric decoherence factor 
| exp(-r plez (t'))\ 2 for an exciton at 77K (bottom panel), as a func- 
tion of A rc , in the presence of a sequence of DD pulses separated by 
At = 0.3 ps for different evolution times: a) t' = 3 ps, b) t' = 30 ps, 
and c) t' = 300 ps. 



i.e. physically, the control cycle time T c becomes significantly 
smaller than the correlation time scale r c = 2ir/u>J£. 

For control separations At = 0.3 ps, u res — 1.04xl0 13 
rads -1 (Tiw res = 6.85 meV), which corresponds to frequen- 
cies uo > lo c such that I exc (id) may become negligible depend- 
ing on the shape of the dot (see Figs. |2]and[3]l. In particular, 
the optimally shaped QD for DD is the one with the small- 
est I exc (uj res ). For V,p e — 455.28 nm 3 and At = 0.3 ps, 
this yields \ rc = 4.8 nm. The top panel of Fig. [9] de- 
picts the decoherence factor due to deformation coupling, 
| exp(— T de f (t'))\ 2 , for an exciton at 4K subject to a DD se- 
quence with At — 0.3 ps evaluated at various evolution times 
t'. Full dephasing may occur for very oblate/prolate QDs, 
however there is also a local maximum in the efficiency of 
the control sequence, which for large t' occurs for dots with 
A ro = 4.8 nm. This is due to this particular shape of dot hav- 
ing the smallest lf x [(id res ). The bottom panel of Fig. |9]shows 
the piezoelectric contribution, | exp(— T piez (t'))\ 2 , for an ex- 
citon at 77K in the same control settings. Once again, the local 
maximum in the DD efficiency for large t 1 is determined by 
the \ rc which minimizes /^f(w res ). At lower temperatures, 
the overall magnitude of the dephasing is reduced. 

The separate components of the spectral density I e , Ih, and 
—2I e h for a QD with A re =4.8 nm (corresponding to the min- 
imum dephasing in the top panel of Fig. [9} are plotted in Fig. 
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[TUl The components of the spectral density include contribu- 
tions from both piezoelectric and deformation potential cou- 
pling. For this reason, the interference term I e h includes both 
negative and positive terms. It can be seen that I exc « Ih for 
uj ^ 8xl0 12 rads -1 . Therefore, close to the local minimum of 
the dephasing we can identify the dot shape needed to mini- 
mize the dephasing by choosing A rc such that Ih is minimized 
at high id. This occurs when the hole wavefunction is spheri- 
cal. For = 455.28 nm 3 and GaAs (i.e. \ rh = 0.495A r J, 
this corresponds to A re = 4.856 nm, which is consistent with 
the minimum dephasing value of A rc = 4.8 nm. The small 
difference is due to the contributions to I exc (oj res ) from I e h 
and I e , see inset of Fig. [10] 

Note that, as expected, as the limit of very fast control is 
approached (T c <C r c ), the shape of the dot has a smaller 
influence on the exciton coherence, which generally remains 
close to its value at the instant the first pulse occurred for a 
long period of time. 



C. DD performance for exciton qubits 

Let us now specifically focus on QDs intended for quantum 
computation applications 6 . We will consider two sets of pa- 
rameters, QD A and QD B. QD A is a GaAs QD with an AlAs 
barrier. The confining potentials are modeled as parabolic in 
all three dimensions, with electron (hole) confinement ener- 
gies in the z-direction of hui e — 505 meV (fkJh, — 100 meV), 
and huj e — 30 meV (Tiidh — 24 meV) in the in-plane di- 
rections, respectively. This corresponds to A re = 6.16 nm, 
\ rh =3.05 nm, and A 2c = X Zh =1.5 nm. QD B is the same 
dot, except an external electric field of 75 kV/cm is applied 
in the x-direction. QD B has been shown to be suitable for 
quantum information processing schemes which exploit direct 
Coulomb coupling between excitons 6 , whereas QD A may be 
useful in schemes not based on the biexcitonic shift. The time 
evolution of the off-diagonal density matrix elements for the 
free dephasing dynamics (Eq. (O) for QDs A and B is de- 
picted in Fig. [TT]for different temperatures 9 . 
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FIG. 11: Free dephasing dynamics of QD A (upper panel), and QD 
B (lower panel) for different temperatures. 



Figure Q~2] shows the qubit coherence evolution for both 
QDs in the presence of periodic DD, according to Eq. ([3]). The 
left column corresponds to QD A, and the right to QD B. Tem- 
peratures are 4K, 77K, and 300K from top to bottom. It can 
be seen that for pulses separated by At < 0.2 ps, DD approx- 
imately freezes the dephasing over the time scales relevant to 
the problem. For a control sequence with longer pulse delays, 
the qubit coherence decays more rapidly, however dephasing 
can still be strongly suppressed for relatively long times. 

In terms of practical quantum computation, the bit-flips do 
not change the time taken for the onset of decoherence, that 
is, the initial time taken for e~ r ^ to sharply decrease from 
unity. The qubit follows its free evolution until the first bit- 
flip occurs, which due to the physical constraint on At, hap- 
pens after some coherence is lost. Still, the pulse sequence 
succeeds at significantly enhancing the exciton coherence at 
large t. For certain operating temperatures and/or QD devices, 
the resulting coherence level may remain very close to unity, 
e.g. for QD B at 4K and At = 0.1 ps, p i(t = 10 ps)= 
V0.997/9oi(t = 0), to be contrasted with the free evolution 
value of %/0.706poi(£ = 0) (see right upper panel). Even at 
room temperature, where for the free evolution all coherence 
is lost in the first picosecond, our calculations show that for 
At = 0.1 ps, poi(i) > \/0.929p i(i = 0) for both QD A and 
QD B for a long period of time. 

It is important to note that in the presence of DD, the qubit 
subject to the applied electric field (QD B) dephases less than 
without applied field (compare left and right column of Fig. 
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FIG. 12: Controlled dephasing dynamics of QD A (left), and QD 
B(right) at 4K (top), 77K (middle), and 300K (bottom), for different 
pulse separations. 



[T2l . This is in contrast to the free evolution behavior, whereby 
the presence of an electric field increases the coupling to the 
piezoelectric field and enhances the dephasing. However, a 
static electric field increases the coupling between excitonic 
qubits and significantly improves the feasibility of quantum 
computation. Our results demonstrate that an electric field 
actually increases the efficiency of periodic DD. 

This improvement may be understood by looking at the ex- 
citon spectral density. Fig. Q~3] compares the spectral density 
for the exciton with and without electric field. It can be seen 
that the spectral density in the presence of an electric field 
is shifted to lower frequencies compared to the no field case. 
This leads to larger long-term dephasing (see Eq. (l26l>). in a 
way similar to the increase of dephasing due to the piezoelec- 
tric field as A re increases. Introducing the electric field ef- 
fectively causes the piezoelectric coupling to contribute more 
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significantly to the dephasing than the deformation potential 
coupling 9 , due to the non-vanishing electron-hole dipole mo- 
ment. Because the relevant spectral density is proportional to 
the Fourier transform of the wave-function, the relative shift 
of the electron and hole wave-functions results in a narrower 
spectral density. In particular, it decreases the spectral density 
at high frequencies (see Fig. [T3T l. As a result, the efficiency of 
sufficiently fast DD sequences (u> c At < ir) is increased due to 
a lower contribution at ui res , as discussed in Sec. III.B. More 
generally, any phenomenon which decreases the spectral den- 
sity at cu res will automatically improve DD performance. 

Also note that for very short pulse separations (u) c At <C ir), 
coherence is maintained near its value at the instant of the 
first control pulse. Since the applied electric field leads to a 
slower initial dephasing than the no field case, less coherence 
has been lost at the instant of the first control pulse, thus poi 
remains at a higher value. Again, this is similar to the piezo- 
electric dephasing behavior observed for increasing X r<s . 

Effectively, the presence of an electric field or a QD shape 
such that X rc is close to X rc , Xf. lez leads to greater dephasing, 
both in terms of the speed of the initial dephasing, as well as 
the magnitude at which coherence saturates, as seen in Sec. 
IIC. However, our results show that it is precisely for these 
QDs that the DD scheme is most efficient, for the shortest 
control sequences allowed by the physical system. Specifi- 
cally, the QDs described in this section have X r<l =6.16 nm, 
which is close to that for which DD is most efficient (see Fig. 
[9). Slightly changing the dot dimensions so that X r<s becomes 
closer to the ideal value would further improve the results re- 
ported in Fig. [12] 



D. DD for short-term quantum storage 

Our calculations indicate that DD may effectively decouple 
the exciton qubit from the environment for a long time, allow- 
ing robust storage of quantum information. In standard solid- 
state quantum computation schemes, holes in double dots 26 
and spin qubits 27 are employed for long-term quantum in- 
formation storage, due to longer T\ lifetimes, whilst exciton 
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FIG. 14: A simple Toffoli gate designed from six controlled-not 
gates, six 7r/8-gates, two Hadamard, and a phase gate. 



qubits are usually proposed for computation only, often as an- 
cillary states, and used for gating purposes. However, quan- 
tum algorithms are generally built up of a series of gates, and 
arbitrary qubit states may have to be stored temporarily while 
other gate operations are carried out. We propose that DD be 
invoked to achieve robust short-term storage. 

As a concrete example, consider the Toffoli gate 1 given in 
Fig. [14] While this is a very simple three-qubit operation, it 
still requires the first control qubit to be stored for the time 
taken to perform 8 gate operations. If we consider an aver- 
age gating time of ~ 1 ps per gate, Fig. [TTI shows that unless 
techniques like the one we propose are used, substantial deco- 
herence would occur on these timescales. 



IV. CONCLUSIONS 

We have shown that the dephasing of an exciton qubit due 
to its coupling with the lattice phonons can be significantly de- 
creased by coherent dynamical control. We found that whilst 
quantum dots with a size and shape most suitable for quan- 
tum computation are the most prone to dephasing, periodic 
dynamical decoupling is most efficient for exactly this type 
of dot. We have demonstrated that a good estimate for the 
shape of the dot which optimizes decoupling performance is 
the one where the hole wavefunction is spherical. In addition, 
our analysis indicates that, although the presence of an elec- 
tric field (which is required by most exciton-based quantum 
computing schemes) increases dephasing, it also makes for 
improved decoupling performance. In the presence of peri- 
odic decoupling, the decay of qubit coherence is substantially 
reduced over evolution times which are very long compared 
to typical gating times. Under favorable circumstances, it is 
reduced so significantly that in practical terms the qubit no 
longer dephases. We propose that dynamical decoupling in 
the simplest periodic implementation examined here may be 
useful to improve temporary data storage during the memory 
stage of quantum computation. 

Besides assessing the impact of pulse imperfections which 
are unavoidably present in real control systems, several ex- 
tensions of the present analysis may be worth considering. In 
particular, dynamical decoupling schemes which involve non- 
uniform time delays and are based on either recursive con- 
catenated design 14 or direct cancellation of high-order error 
terms 28 have been recently found to attain remarkable high 
decoupling fidelity under appropriate assumptions 14 ' 18 ' 29 . Fur- 
thermore, so-called Eulerian decoupling schemes 30 allow de- 
coherence suppression to be achieved without requiring unre- 
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alistic control strengths as in the bang-bang limit - which may 
be critical to minimize unwanted excitations. A future investi- 
gation will be needed to identify the possible added benefits of 
more elaborated decoupling schemes under the specific phys- 
ical constraints and design trade-offs associated with exciton 
based quantum-dot devices. 
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